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Abstract 

Applying a technique developed in a recent work[l] to calculate wavefunction evolution in 
a dissipative system with Ohmic friction, we show that the wavelength of the wavefunction 
decays exponentially, while the Brownian motion width gradually increases. In an interference 
experiment, when these two parameters become equal, the Brownian motion erases the fringes, 
the system thus approaches classical limit. We show that the wavelength decay is an observable 
phenomenon. 
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1. Introduction. The simplest example of a dissipative system, an harmonic oscillator 
coupled to an environment of a bath of harmonic oscillators, has been the subject of extensive 
studies [2-4] (see [1] for more references). In a recent paper [1], we obtained a simple and exact 
solution for the wave function of the system plus the bath, in the Ohmic case. It is described 
by the direct product in two independent Hilbert spaces. One of them is described by an 
effective Hamiltonian, the other represents the effect of the bath, i.e., the Brownian motion, 
thus clarifying the meaning of the wave function of the effective Hamiltonian after the bath 
variables are eliminated, establishing a connection between two different approaches to the 
dissipative system problems: the effective Hamiltonian approach, and the system plus bath 
approach. 

In this paper we shall apply this result to a specific case to show how the wavelength of the 
wavefunction decays exponentially while the Brownian motion width increases at the same 
time. We shall show that in an experiment with two overlapping wave packets interfering 
with each other, the interference term of the probability density decreases rapidly when the 
wavelength of the wavefunction of the effective Hamiltonian reduces to be equal to the Brownian 
motion width. Therefore, our result provides a transparent picture about how a quantum 
system, when interacts with environment and hence becomes a dissipative system, gradually 
evolves into a classical system. 

First, we derive a simple formula for the probability density. Next, we apply our method to 
the case of an initial gaussian wave packet displaced from the center of the parabolic potential 
well, and show that the wavelength of the effective Hamiltonian decreases exponentially while 
the Brownian motion width increases. Then, we use the results to calculate the interference 
term, assuming the initial state is a superposition of two wave packets with minimum un- 
certainty, one of them is centered at the origin, the other is displaced from the center. This 
problem has been studied by A.O.Caldeira and A.J.Leggett (CL)[4]. However, the result does 



2 



not show the expected fringe wavelength decay. This puzzle stimulated our further analysis 
given in the last section, which explains why the wavelength decay is not observable in this 
well studied case, and shows that if the initial wave packet width is different from that with the 
minimum uncertainty, the wavelength decay becomes an observable phenomenon . We shall 
use the same notation as the previous paper, and quote formulas from there. 

2. The probability density. We consider an harmonic oscillator of frequency ujq, linearly 
coupled to a bath of harmonic oscillators of frequency cuj with Ohmic spectral density distri- 
bution. The Hamiltonian of this system is given in reference [1] : 

It is well known that the coupling to the bath introduces damping to the harmonic oscillator 
[1]. The damping rate is rj and the system frequency is shifted to a; = {lOq — rf /Ay/'^ by the 
damping. In the previous paper [1] we showed that the wavefunction of the system plus the 
bath, in the Schoedinger representation, at time t, can be written as: 

N 

*(g,{0},i) = V'(g-EO'^) JlUij.t) , (2) 

i i=i 

where the wavefunction ip{Q,t) is a solution of the effective Hamiltonian eq.(23) of [1]: 

= + ^Afa;o^e''*Q' , (3) 

while — bji(t)xjo + bj2(t)xjo is the contribution of the j'th bath oscillator to the Brownian 
motion with Xjo,XjQ its initial position and velocity. The coefficients bji{t) and bj2{t) are well 
known in elementary physics[l]. The function XjiCj^'t) is given by XjiCj^'t) —< %|XjO 
where Xjoi^jo) is the initial state of the bath oscillator, and 16*^^ > is an eigenstate of the 
operator with eigenvalue ^j, determined to within an arbitrary phase, which depends on 
only ^f. 
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(4) 



We are interested in the probability density time evolution if initially the bath is in an 
equilibrium at temperature T. Assuming initially the j'th bath oscillator is in an excited state 



Uj, we write the initial wavefunction of the this oscillator as Xj i^jo) 



X 



(nj) 



Xjo). Then, the 



probability density for finding the system at position q and finding the contribution of the 
bath oscillator j to the Brownian motion to be ^j, is: 



N 



with the distribution of at time t: 



(5) 



where /? = i/kT, and the partition function is 



(6) 
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n,=0 



2sinh(^) 



(7) 



Using Xji^jy i) =< % \Xjo >, we find Pj{Cj] ^) =< % \Pjo\ % where pjo is the initial density 
matrix of the j'th bath oscillator. The matrix elements of pjo is well known in quantum 
statistics text book[3]: 



Pjoi^jO^^'jo) — 



rrijUJj tanh( 



phujj 



Exp 



- coth.((3hoi)j){x^Q + 



XjoX 



smh.(P?iujj 



(8) 
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Hence, using eq.(^ , followed by a straight forward calculation, we obtain (obviously the 
arbitrary phase of 9^_. would not influence the result): 
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which is a gaussian distribution with width: 



a 



(|6,i(t)|^ + a;;|6,2(t)|^)coth — ^ . (10) 



The density matrix traced over the variables is given by: 

N 1 _ -J 

p(gi,g2;t) = / / /.../r(gi-Eo,^M?2-Eo,^)n7|— e'^de,. (ii; 



N 

Notice that this is not a reduced density matrix in the ordinary sense, because the meaning 
of its off diagonal elements is not clear, mainly due to the fact that the variable ^j(t)'s are not 
pure environment variables unless t=0. However, the diagonal matrix elements p(g, g; t) have 
clear physical meaning: they are the probability density p{q;t). Hence we shall only address 
the diagonal elements in the following, even though the results for the off-diagonal elements 
are as simple as the diagonal elements. 

It is easy to derive the following formula for an arbitrary function f{C,)'- 

///...//Ean "^"(A-^Jme (12) 



N 

with 



N 

4 = E^l- (13) 



Therefore we have a very simple expression for the probability density: 



p(g,t) = ^=^ / |^(g-e,t)re"^de (14) 
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where is the Brownian motion width, calculated by substituting eq.(p!OD into eq.( |T3|): 

alit) = I ^^i\b,^it)\' + co]\b,,{t)f) cothi^\p{u;,)du;,, (15) 

and by replacing the sum over the bath oscillator number j with the spectral density p{ujj) of 
the bath oscillators as given by eq.(6) of [1]: 



2r]MmjUj] 

Pi^j) = — g^- (16) 

This width is zero initially, then approaches its final equilibrium value in a time interval 
of the order of l/rj. In low temperature limit, the equilibrium width is simplified to: 

h 



a^(t = oo) 



2nMLJ 



— h arctan — 
2 \r]uj^ 



(17) 



If <^ ujq, this width happens to become the same as the width of the ground state of the 
system cTq = %/ {2Mujq). 

There are some subtleties to be addressed in future communications in the calculation of 
cr^(t) for finite t. They are associated with a logarithmic divergence of the integration over 
ujj in eq.(15). The divergence is removed by introducing a cut-off frequency cuj'^"*'' in eg . (|T5|) . 



In the following, for simplicity, we consider only the low temperature limit. To get an idea of 

t|/o"q in figure 1, for uj^^ 



the behavior of a^, we plot o'ljo'l in figure 1, for tuj'^"*^ = lOOu;, and r/ = O.lcuo- As a rough 



estimate, we have 



This function is not very sensitive to the specific value of cuj^^^/cjo- 

3. The wavefunction evolution. We assume that initially the wavefunction is a gaussian 
wave packet with minimum uncertainty and displaced from the center of the potential well by 



^o(g)=^(g;t = 0) = (27ra2)-i/^e _ (^g) 



Using the Green's function (derived as eq.(21) in [1] ): 



G{Qi, Qo; t,0)= I — — I exp 



iM 
2ha2 



we have: 



i^{q,t)= / i'{qo,0)G{q,qo;t,0)dqo 



(27r<7o^)-V^ 



Voi + iuJoa2 \ 4(7o(ai + iuj^a^) 
where Oi and are the well known coefficients of the initial values and for the solution 

of the classical damping equation q + rjq + uj^q = [1]. If 77 -C uJqi we rewrite ai and 02 

(eq.(12) of [1]) approximately as: 



0,2 + iuoa2)e^'*q^ — 2qz + aiz^ 



(20) 

(21) 

(22) 



V 



ai — e 2*(cosa;t + 7^ sina;^) fa e ^^coscut 
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uJoa2 — e 2 — sinut^e sin cot, d2{t) ^ ai{t). 

OJ 

With this approximation, the wave function takes a simple and familiar form: 



(23) 
(24) 



ilj{q,t) ^ C{t)e 4(72e ^iK{t)q_ 



(25) 



This is a gaussian wave packet centered at the classical damped orbit: qc{t) — e z cos ujt, 
with momentum of 

hkJt) ^h-^e^^sinutc^ Me"^ ^qJt) . (26) 
2(To [at J 

This corresponds to the wave number of a particle of mass Me***. The coefficient C{t) is 
simply a normalization constant (independent of q). Notice that even though the velocity 
of the classical motion qc{t) deceases as e"''*/^, the wavelength 27r /kc{t) still decays as e"''*/^ 
because the effective mass increases exponentially. Notice also that the width of the wave 
packet (Toe"''*/^ deceases exponentially, so that the probabihty density distribution evolves 
into a 5 function if Brownian motion can be neglected. Obviously, as time approaches infinity. 



the system approaches a classical limit. To provide a more clear picture of the evolution 
process, in Figure 2, we plot the real part of the wavefunction at 4 different times, assuming 
that the damping rate is 1] = O.Iuq, and the initial displacement is z = Sctq. 

4. The interference fringe wavelength. Now, we are ready to check if the wavelength 
decay is observable. If initially the wavefunction is a superposition of two such wave packets 
as eq.( 1T9|) , one with z=0, the other with zj^O, the above results can be used to calculate 
straightforwardly the interference term (which has been studied in detail by CL using path 
integral technique [2,4]). The result will show how the Brownian motion erases the interference. 

The initial wavefunction is: 



where C is a normalization constant. 

Using eq. (p7D , eq. ([2^) , and eq. (p^ , we find: 



p(g,t)=C7^- 



e'^e + e"^^ + e ^2 cos ( l^lii^^^^ ] ^"^^ 



f28) 



(ye 
where 

al = al{a\ + ujlal)+al (29) 

The last term in the parenthesis is the interference term, which agrees with the equation (2.9) 
of [4]. The notation is slightly different. For example, a part of the exponent of the last factor 
has been identified as the ratio of the Brownian motion width cr| over the full width cr^ of 
the wave packet; ai is identified as the trajectory of a particle with initial position 1 but zero 
velocity, while 02 corresponds to a particle with initial velocity 1 but position at the origin [1]. 
In every period, the maximum interference occurs when ai = 0, {ut ^ {2n + l)7r/2), as is 
indicated by the first factor of the interference term. In the following, we only discuss the 
time when ai(t) = 0. 
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Since the wavelength decays exponentially in eg. p^) , it is interesting to see whether the 
interference fringe wavelength also decays exponentially, and provides an observable phe- 
nomenon. We assume r] <^ uq, use the approximation eg. (p4D , and estimate the Brownian 
width by eg . ([ISD . Then the growth of o"| and the decay of happen to cancel each other, 
and (jg ^ ctq is nearly constant. The eg.(^) becomes: 



_ ,2 

p{q, t) ^ 2C^e [l + cos (fccoe"^^*?) e"^'^*'"!] , (30) 

where fcco = ~ kc{ujt = ^) (see eg.(^U])). In order to have enough fringes the initial 
wavelength 27r//cco should be much smaller than the wave packet width a^. Hence the Brownian 
width increases to the initial wavelength in a time much shorter than the damping time l/rj 
(the time for it to increase to ctq ). Therefore the eg. (|30|) shows that the fringe wavelength does 
not decay but increases slightly before the Brownian motion erases the interference pattern. 
Thus, in this well studied case, the exponential decay of wavelength is not experimentally 
observable! 

5. How to observe the wavelength decay. If we carefully examine the calculation of the 
probability density eg . (|I^) , we realize that the slight increase of fringe wavelength is due to 
the narrow wave packet width ctq and the convolution of the wavefunction with the Brownian 
motion in eg.(0). To see whether the wavelength decay is really observable, in the initial 
wavefunction eg.(^) we replace cxo by a different width a. Let the ratio <jQ/a = r. Now, with 
a straightforward calculation as before, the probability density becomes: 



+ e~~^ + e ^2 cos ( S^zllJZ^^A e"^^ 



while 



(31) 



a^g = a\al + rW^al) + (t|. (32) 



Again, we use eg.(|T8D to calculate the fringe wavenumber (the coefficient of g in the phase of 



the cosine function in eg. (|3T|) ) . The fringe wavenumber is then found to grow exponentially 
until t ^ /n(r^ — l)/ri, before it starts to decrease, due to the Brownian motion and the finite 
wave packet width. In the previously studied case, r=l, hence the wavelength decay is not 
observable. However, if r is sufficiently large, it becomes observable, as can be seen as follows. 
Now, for quite a long time, ct| can be neglected in eq. (|32D , while ag decays exponentially, 
the interference pattern can then be simplified to (when ai = and when the approximation 
eq.dl) is valid): 



2a 



q 



p{q,t) = C''—e 1 + cos e-^'^'^^H . (33) 

In the mean time, the visibility of the fringes also decays. The time when it is reduced 
by a factor of e is determined by -^-^ — - = 1, i.e., when the wavelength is about equal to 
the Brownian motion width. Using eq. ([I8|) , we find the time is t ~ ln{l + Sctq'^/z'^)/')] . 
Thus large r and small z tend to give longer time to observe the wavelength decay before the 
fringes disappear. The number of visible fringes in the wave packet is 2kc{t)ag 'r^ zja when 
t {2n + l)vr/2, so z should be larger than a to have enough fringes. In figure 3 we plot 
the probability density at different times for a case with z = 3aQ and r =16. A comparison 
of fig. 3b and fig. 3c clearly shows the wavelength decay and the visibility decay. In fig.3d, the 
coherence has been erased almost completely, long before the distribution width becomes ctq. 
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Figure Captions 



Figure 1. The Brownian motion width as a function of time. The cut-off frequency 
Figure 2. The time evolution of ip{q,t). z = Sctq,?] = O.Iuq. 

Figure 3. The time evolution of the interference pattern. The probability density is in 
arbitrary unit, z — Sctq, rj — O.luo, r — 16, o;]™*^ = lOOo;. 
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